Molecular CO(l-O) gas in the z ~ 2 radio galaxy MRC 0152-209 
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ABSTRACT 

We report the detection of molecular CO (1-0) gas in the high-z radio galaxy MRC 0152-209 
(z = 1.92) with the Australia Telescope Compact Array Broadband Backend (ATCA/CABB). 
This is the third known detection of CO(1-0) in a high-z radio galaxy to date. CO(1-0) is the 
most robust tracer of the overall molecular gas content (including the wide-spread, low-density 
and subthermally excited component), hence observations of CO (1-0) are crucial for studying 
galaxy evolution in the Early Universe. We derive L' co = 6.6 ± 2.0 x 10 10 K km s _1 pc 2 for 
MRC 0152-209, which is comparable to that derived from CO(1-0) observations of several high-z 
submillimcter and starforming BzK galaxies. The CO(1-0) traces a total molecular hydrogen 
mass of M H2 = 5 x 10 10 (a x /0.8) M . MRC 0152-209 is an infra-red bright radio galaxy in 
which a large reservoir of cold molecular gas has not (yet) been depleted by star formation or 
radio source feedback. Its compact radio source is reliably detected at 40 GHz and has a steep 
spectral index of a — —1.3 between 1.4 and 40 GHz (4 — 115 GHz in the galaxy's rest-frame). 
MRC 0152-209 is part of an ongoing systematic ATCA/CABB survey of CO(1-0) in high-z radio 
galaxies between 1.7 < z < 3. 

Subject headings: Galaxies: active — Galaxies: evolution — Galaxies: high-redshift — Galaxies: indi- 
vidual (MRC 0152-209) — Galaxies: ISM — Radio lines: ISM 

1. Introduction 



High-redshift radio galaxies (HzRGs; Lsoomhz > 
10 27 W Hz^ 1 ) are among the most massive galax- 
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iesin the early Universe (jMilev fc De Breuck 
2008). They are typically surrounded by rich 
proto-cluster environments and are believed to be 
the ancestors of early-type central-cluster galaxies 
fe.g. IVenemans et al]|2007l ). HzRGs are in an ac- 
tive stage of their evolution, with powerful radio 
jets emanating from a massive central black- hole. 
These radio jets vigorously int eract with their 



surro unding environment (e.g. iHumphrey et al 



2006) and serve as beacons for tracing the faint 
host galaxy and surrounding proto-cluster. This 
makes HzRGs among the best studied high-z ob- 
jects and ideal laboratories for investigating both 
the formation and evolution of galaxies/clusters 
as well as the relationship between early star for- 
mation and AGN activity. 

A crucial component in the evolutionary studies 
of HzRGs and their associated star formation pro- 
cesses is a detailed knowledge about the content 
and properties of the molecular gas. The most 
abundant molecule, molecular hydrogen or H2, 
has strongly forbidden rotational transitions and 
is virtually invisible, except when it is shocked or 
heated to high temperatures. An excellent tracer 
for H2 gas is carbon monoxide or CO, which emits 
strong rotational transition lines that occur pri- 
marily through collisio ns with H2, even at rela- 



tively lo w densities (e.g. ISolomon fc Vanden Bout 
20051 ). iBrown fc Vanden Boutl (|l99ll ) were the 



first to observe 12 CO (referred to as CO in this 
paper) beyond z = 2. 

During that same decade, intensive searches 
failed to detect CO in HzRGs (|Evans et al.lll996l 

Since then, CO emis- 
individual studies 
2-5 



van Oiik et all |l997a) 



sion has been detected in 
of powerful radi o galaxies between z 
see reviews bv ISolomon fc Vanden Boutl 120051: 
Milev fc De Breuckl 120081 and Sect. I4T21 for in- 
dividual references). In some cases the CO has 
been resolved on scales of several tens of kpc (e.g. 
Papadopoulos et al. I l2000h . 

Despite these results, two major limitations to 
systematic searches for CO in HzRGs have been 
the very limited velocity coverage of existing mm- 
spectrometers (often not much wider than the ve- 
locity range of the CO gas and/or the accuracy of 
the redshift) and the fact that most observatories 
can only target the higher rotational transitions of 
CO at high-z. While these higher CO transitions 
trace the dense and thermally excited gas (such 



as that in the nuc le ar st a r burst /AGN region), 
Papadopoulos et al.l (|2000l . l200ll ) suggest that 
more widely distributed reservoirs of less dense 
and sub-thermally excited gas could be largely 
missed by these observations and could thus more 
easily be detected in the low er transitions than 



often assumed (se e also e.g. |D anncrb auer et al . 
20091 : ICarilli et al.l l2010l : llvison et al.l l201ll and 



references therein). The luminosity of the ground 
transition, CO(1-0), is least affected by the excita- 
tion conditions of the gas. Therefore, observations 
of CO (1-0) are crucial for deriving the most accu- 
rate estimates of the overall molecular gas content 
in high-z galaxies. 

Ongoing developments in broadband correla- 
tors are overcoming the ab ove mentioned limita- 
tions. As we discussed in lEmonts et al.l (|201ll ). 



the upgrade of the Australia Telescope Compact 
Array (ATC A) with the 2 x 2 GHz Compact Array 
Broadband Backend (CABB; IWilson et ai1l201lh 
now allows searching for CO (1-0) across the south- 
ern hemisphere. We are currently in the process of 
observing an unbiased sample of HzRGs from the 



Molonglo Reference Catalogue ([McCarthy et al 



1990h in CO(1-0) with the ATCA/CABB 7mm 
system between 1.7 < z < 3. 

In this Letter, we present the detection of 
CO(1-0) emission in one of our sample sources, 
MRC 0152-209 (z=1.92). MRC 0152-209 is a ra- 
dio galaxy with a dis turbed optical morphology 
( Pentericci et al.l 200lh. It is bright in the near 



and m id-IR ( Seymour et al.l2007t De Breuck et al 



2010h . detected in Lva ((McCarthy et al.l Il99ll ) 



and has a co mpact radio source (1. 6 arcsec/13 kpc 



in diameter; Pentericci et al. 2000h 



Throughout this paper, we assume Ho = 
71kms- 1 Mpc~ 1 , n M = 0.3 and fi A = 0.7, corre- 
sponding to an angular distance scale of 8.3 kpc 
arcsec -1 and luminosity distance Di, = 145 71 
Mpc for MRC 0152-209 (following [Wrigl3[2006) 



2. Observations 

MRC 0152-209 was observed on 25 & 26 August 
and 28 & 29 September 2010 with ATCA/CABB 
in the most compact hybrid H168 and H75 array 
configurations. The total effective on-source inte- 
gration time was 15. 5h (with 45% of our overall 



1 See http://www.astro.ucla.edu/~wright/CosmoCalc.html 



2 



observing time spent on overheads and calibra- 
tion). We used the coarsest spectral- line CABB 
mode (2x2 GHz bandwidth; 1 MHz resolution) 
and centered both 2 GHz observing bands around 
39.476 GH^], the frequency that corresponds to 
the CO(l-O) line (V rcst = 115.2712 GHz) Doppler- 
shiftcd to z — 1.92. We observed only above 
30° elevation and during atmospheric conditions 
in which there was no significant phase decorre- 



lation of the signal, with T s 



50 - 80K dur- 



ing the various runs. Short (2 min) scans were 
made on the nearby strong calibrator PKS B0130- 
171 (54o GHz ~ 2.3 Jy beam -1 ) roughly every 
10 minutes in order to derive frequent gain- and 
bandpass-calibration solutions, both of which were 
interpolated in time and applied to the data of 
MRC 0152-209. We observed Uranus when it was 
at roughly the same elevation as PKS B0130-171 
in order to flux calibrate our data with an accu- 
racy of 30%. Details on the calibration strategy 
(using miriad) and flux calibration uncertainties 
are described in lEmonts et al. ( 2011 ). 

We used multi-frequency synthesis to o btain a 
conti nuum image (with robustness = +1; iBriggsl 
19951 ). which we subsequently cleaned. We then 



separated the spectral-line from the continuum 
data with a linear fit to the line-free channels 
in the wy-domain before Fourier transforming the 
data into a line-data cube (robustness = +1). 
The velocity axis was translated to match the op- 
tical barycentric rest-frame vel ocity at the red- 



shift of MRC 0152-209 (z = 1.92; lDe Breuck et al 
2000). The spatial resolution of our final data 
is 9.80 x 7.09 arcsec (PA = 89.8°), with a veloc- 
ity resolution in the line data of 7.6 kms -1 per 
1 MHz channel at 39.476 GHz across an effec- 
tive velocity range of ^15000 kms -1 . The noise 
is 58 /zJy beam -1 in the continuum image and 
0.47 mjy beam -1 per 1MHz channel in the line 
data. For the data analysis presented in this pa- 
per we have binned the line data to 40 kms -1 
wide channels and subsequently applied a han- 
ning smooth, resulting in an effective velocity res- 
olution of 80 kms -1 and a noise level of tr = 0.16 
mjy beam -1 chan -1 . 



2 The actual central frequency of the bands was chosen to be 
39.9 GHz, because channels were missing in the first half 
of each band at the time of our observations. Only one of 
the observing bands was used in the final data analysis (sec 
also lEmonts et"aDl201 if ) . 
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Fig. 1.— Top: CO(1-0) spectrum of MRC 0152- 
209 against the central host galaxy (extracted by 
averaging the central four pixels of the binned line 
data, Sect. [2]). Bottom: zoom-in of the CO(1-0) 
spectrum. The dotted line show a Gaussian fit. 

3. Results 
3.1. CO(l-O) 

Figure Q] shows the spectrum of CO(1-0) against 
the central region of MRC 0152-209. From fit- 
ting a Gaussian to this profile, we derive that the 
CO (1-0) signal has a peak flux density of Sco = 
0.84 ± 0.25 mjy beam -1 E and FWHM = 397 ± 37 
kms -1 . The CO(1-0) signal peaks at a redshift 
of zco = 1.9212 ± 0.0002, in agreement with 
the optical redshift of z = 1.92 ([McCarthy et al. 
199ltbe Breuck et al.ll2000l ). FigureE] (left) shows 



d The error reflects the accuracy of the flux calibration (Sect. 
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Fig. 2. — Left: total intensity image of the CO(l-O) signal of Fig. [T] (see text for details). Contour levels: 
-3, -2 (white); 2, 3, 4, 5, 6, 7 (black) a, corresponding to L' co = -2.5, -1.7 (white); 1.7, 2.5, 3.3, 4.2, 5.0, 
5.8 xlO 10 K km s _1 pc 2 . The inset shows the same CO contours overl aid onto the unresolve d 40 GHz 
ATCA radio continuum (grey-scale). Right: HST/NICMOS image from IPentericci et all d200lh . showing 
the di sturbed host-galaxy morphology. The inset show the 8.2 GHz radio continuum from IPentericci et al.l 
(|2000h . 



a total intensity map of the CO(l-O) line inte- 
grated over the channels in which the line was 
detected. The integrated CO(l-O) signal is de- 
tected at the 8cr level and / v S C oSv = 0.37 ± 0.11 
Jy beam -1 x krns" 1 . There is a tentative indica- 
tion from Figure [5] (left) that the CO(1-0) emis- 
sion is marginally resolved in west-north-western 
(WNW) direction. This would mean that the 
CO(1-0) gas is distributed on scales larger than 
the major beam axis, which corresponds to ~81 
kpc. Observations at higher spatial resolution are 
in progress to verify this. 

The CO(1-0) emiss ion-line luminosity can 
be calculated following Solo mon fc Vanden Boutl 
([20051 . and references therein): 



L'co = 3.25xl0 7 (^^)(^) 2 (^)- 2 (l + ,)- 1 , 
Jy km/s Mpc GHz 

with L' co expressed in K km s 1 pc 2 . For MRC 0152- 
209, L' co = 6.6 ± 2.0 x 10 10 K km s" 1 pc 2 . 

3.2. Radio continuum 

The radio continuum emission of MRC 0152- 
209 is unresolved in our ATCA observations and 
has a flux of S39. 9GHz = 5.1 ± 1.5 mJy Using our 
40 GHz flux and 8.2, 4.7 and 1.4 GHz data from 



Pentericci et al. (2000), the spectral index can be 
fitted with a power-law of the form S v oc v a from 
~1.4 to 40 GHz (4-115 GHz in the galaxy's rest- 
frame; Fig. This yields a steep high-frequency 
spectral index of a4°4GHz ~ — 1-3- Data obtained 
from surveys conducted below 1 GHz suggest a 
flattening of the spectrum at low frequencies, pos- 
sibly fitted more accurately by a second order 
polynomial rather than a single power-law (Fig. 
|3|). Additional measurements are required to de- 
termine this reliably, but this is beyond the scope 
of the current paper. 

4. Discussion 

4.1. Molecular gas and star formation 

From the CO(1-0) luminosity we can esti- 
mate the mass of molecular gas in MRC 0152- 
209 by adopting a standard conversion fac- 
tor a x = M H2 /L^ [M© (K krns -1 pc 2 )" 1 ] 
(where Mh2 includes a heliu m fraction; e.g. 
Solomon fc Vanden Boutl 120051 ). In this paper, 
we use a x = 0.8 (Kkms -1 pc 2 ) -1 found 
by iDownes fc Solomonl ( 19981) for ultra-luminous 
infra-red galaxies (ULIRGs). This also agrees 
with a x = 0.8 - 1.6 M (K kins" 1 pc 2 )" 1 found 
for high-z sub-millimeter and star forming galax- 
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3.— Radio spectral index of MRC 0152 
Data at 1.4 4.7 and 8.2 GHz 
(l2000h 



from Pcntericci et al 



are 

the lower frequency 



data are taken from various surveys listed in 
the NASA Extra-galactic Database (NED - see 
|http://nedwww.ipac.caltech.edu/[ ). The solid line 
represents a broken power-law fit {F v oc v a ), the 
dashed line shows a second order polynomial fit. 



ies (|Tacconi et al.l 120081: IStark et ail 120081) . We 

note, however, that the conversion of CO lumi- 
nosity into H2 mass depends on the properties 
of th e gas, such as metallicity and radiation field 



(e.g. 



Glover fe Mac Low! 120111) . and that conver- 

5 Mr (K kms- 1 



sion factors as high as of a x = u iviq 
pc 2 ) -1 have been deri ved for molecular clouds i n 
the Milky Way (e.g. ISolomon fc Barrettl fl99ll ). 
Our assumption of a x = 0.8 M Q (K kms -1 
pc 2 )" 1 thus results in a conservative estimate of 
the total molecular gas mass in MRC 0152-209 of 
M H2 = 5 x 10 10 M Q . 

The 24/xm infra-red (IR) flux of MRC 0152-209 
(/24^tm = 3320 ± 133 /LtJy) is among the high- 
est in a l arge sample of HzRGs observed with 
Spitzer by |Pe Breuck et al.l (|2010f ). At z = 1.92, 
the observed A = 24 (tun corresponds to A = 8.2 
/im in the rest-frame, hence from /24/jm we de- 
rive a mono-chromatic rest-frame luminosity of 
(vLuYlf^ = 2.7 x lO 12 !, ^. When comparing the 
mid-IR rest- frame colors ( De Breuck et al.l 2010l) 
with IRAC color diagrams bv lStern et al.l u2005l) . 
MRC 0152-209 falls within the region that is dom- 
inated by active galaxies, indicating that a signif- 
icant fraction of the 8 /zm flux may originate from 



a hot dusty torus surrounding the central AGN. 

Alternatively, part of {vL v ) T ^ may arise 
from a dust-enshroud ed starburst, which peaks 
in the far-IR (FIR). lEvans et all (|2005l) show 
that low- z radio galaxies follow a strong corre- 
lation between L' co and Lfir- This trend is also 
seen among high- z submillimeter gala xies (SMGs; 



Greve et al.l 120051 : llvison et al.l 120111 ) and likely 



reflects the Schmidth-Kennicutt relation between 
star formation rates and gas reservoir (|Schmidd 
19591: " 



Kennicutt 19981) . If this correlation 



is 



also valid for HzRGs, the CO(1-0) luminosity of 
MRC 0152-209 implies a FIR luminosity of or- 
der L Fm - 5 x 10 12 L Q (llvison et al.ll201ll). wel l 



Bavouzet et al.l 12008) 



in the regime of ULIRGs 
show that the mono-chromatic IR luminosities 
of starburst-dominated ULIRGs at 8, 24, 70 and 
160/J,m are strongly correlated with the total IR 



luminosity (£ir = £5-1000^™): with reliable cor- 
relations up to z ~ 2. While these correlations 
may not necessarily be valid in case a significant 
fraction of (uL v ) T s % m in MRC 0152-209 is AGN 
related, we argue that they still provide a reli- 
able upper limit to Lir, given that a 'starburst- 
only' scenario will provide the highest possible 
energy output in the FIR and thus maximize 
-L tr- F rom {vL^)™^ and following lBavouzet et al 



(|2008l their Equation 6), we therefore estimate 
that L m < 7.9 x 10 12 i Q for MRC 0152-209. 

The upper limit on the dust-to-gas ratio 
in MRC 0152-209 can then be estimated to 
be Lir/L^q < 120, consistent w ith that of 
nearby ULIRGs JSolomon et al.lll997l ). Following 
Kennicuttl ( 19981 ). the upper limit on the star for- 
mation rate is SFR » 5 S xw'->l q ~ 1362 M J* -1 - 
Assuming that the molecular gas in MRC 0152- 
209 offers a reservoir to sustain this maximum star 
formation rate, the lower limit on the gas depletion 
time-scale is t dc pi ction = |^ > • 39 Myr (with 
< e < 1 the sta r formation efficiency). Ac cord- 
ing to models by iMihos fc Hernquisd ( 1994[) . this 
minimum gas depletion time-scale is comparable 
to the typical lifetime of a massive burst of star 
formation during a major merger, which could give 
the merging system the appeara nce of a ULIRG 
for about 



50 Myr (although Swinbank et al 



(|2006l ) show that for high-z SMGs a starburst 
age of several 100 Myr is more realistic). It is 
conceivable that also MRC 0152-209 has under- 
gone a recent merger, given its disturbed optical 
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morphology (jPentericci et al.ll200ll) . 

Therefore, MRC 0152-209 shares characteristics 
commonly found among ULIRGs, but its classifi- 
cation remains somewhat ambiguous until the na- 
ture of the IR emission and star formation prop- 
erties are better understood. It is clear, however, 
that MRC 0152-209 contains a large reservoir of 
molecular gas that has not (yet) been depleted by 
vigorous star formation. 

4.2. CO(l-O) in HzRGs 

The CO(l-O) luminosity of MRC 0152-209 is 
an order of magnitude larger than that found in 
nearby radio gala xies ( Ocana Flaauer et al. 2010t 
Evans et alj|2005l) . indicating that MRC 0152-209 
contains significantly more molecular gas than its 
low-z counterparts. 

Assuming a x = 0.8, the H 2 mass of MRC 0152- 
209 is of the same order as t hat derived fro m CO 
detections in other HzRGs rtScoville et al 



1997 



Papadopoulos et al. 2000, 2001; Alloin et al. 2000; 
De Breuck et alj|2003afbl booa Greve et al.H2004 ; 



IJ\ L»1LUU\ LL 0,1. ZiUUUa, U , ^ W W 'J , VJl LVL dl. 

Klamer et al.l2005Hlvison et al.l2008tlNesvadba et al 



2009 ) as well as f rom a survey of high-z SMGs 



dGreve et al.ll2005h . However, most of these CO 



studies relied on observations of the higher rota- 
tional CO transitions, which could underestimate 
the total mass and spatial extent of the molecular 
gas (see Sect. [1]). 

To our knowledge, CO(l-O) detections have 
been claimed for only t wo other HzRGs, na mely 
TN J0924-2201 (z = 5.2:lKlamer et al.ll2005l) and 
4C 6 0.07 (z = 3.8; iGreve et al J 12004 llvison et all 
20081 ). Both systems contain M H 2 ~ 10 n M Q . 
Interestingly, they both are IR-faint objects 
(|De Breuck et al.ll2010t llvison et al.ll2008l ). CO(l- 
0) observations o f four SMGs (z = 2.2 - 2.5; 
Ivison et al.ll2011 ) and three massive star-formin; 
disk (BzK) galaxies (z ~ 1.5; Aravena et al.ll201 
reveal CO(l-O) luminosities comparable to that of 
MRC 0152-209. CO (1-0) observations of larger 
and unbiased samples of various types of high-z 
galaxies are necessary to investigate and compare 
the molecular gas content in these systems in a 
systematic way. 

4-. 2.1. CO(1-0) and radio source properties 

The radio source in MRC 0152-209 is fairly 
compact (13 kpc) and not reaching beyond the op- 



tical bound aries of the host galaxy (jPentericci et al 
2000L l200lh . Based on a CO(3-2) detection in the 
outer part of the Lya halo that su rrounds the z = 
2.6 ra dio galaxy TXS 0828+193, iNesvadba etal 



(2009) speculate that large amounts of dense gas 
(and perhaps dust) may be present in the halos 
of HzRGs at radii that are not yet affected by 
mechanical heating from the radio source. This 
would be in agre ement with studies of t he neutral 
gas in HzRGs bv lvan Oiik et al.l (|l997bl ). who find 
that strong H I absorption is predominantly asso- 
ciated with small (< 50kpc) radio sources. It will 
be particularly interesting to investigate whether 
the CO (1-0) gas, i.e., the part of the observable 
molecular gas that is least affected by excitation 
and heating, shows a similar trend with radio 
source size, or whether on average HzRGs are 
fainter in CO(l-O) than their radio-quiet high-z 
counterparts. 

Besides MRC 0152-209, we recently published 
a tentative CO (1-0) detectio n with ATCA/CABB 
in the HzRG MRC 0943-242 (jEmonts et al.ll201ll) . 
This tentative CO(1-0) detection is found in 
the outer part of the Lyct halo that surrounds 
MRC 0943-242 (i.e. beyond the -15 kpc ra- 
dius of the radio source). Regarding the other 
two known HzRGs with CO(1-0) detections (Sect. 
H2J, TN J 0924-2201 also cont ains a compact ra- 
dio source (|Klamer et al.ll2005l ). while for 4C 60.07 
the radio source is extended, but the CO (1-0) most 
likely represents tida l debris from a gal axy-galaxy 
interaction (z = 3.8; Ivison et al. 20081) . 

In a forthcoming paper, we will present result of 
our ongoing systematic CO(1-0) survey of a larger 
sample of HzRG with ATCA/CABB. This will al- 
low further investigation into possible connections 
between radio- AGN activity and the cold molecu- 
lar gas properties of HzRGs. 

5. Conclusions 

We have presented a robust detection of CO(l- 
0) in the high-z radio galaxy MRC 0152-209 (z = 
1.92), tracing a molecular gas mass of Mh2 = 5 x 
10 10 (a x /0.8) M . MRC 0152-209 is an IR-bright 
radio galaxy that contains a large reservoir of cold 
molecular gas that has not (yet) been depleted by 
star formation or radio source feedback. 

CO (1-0) is the most robust tracer of molecular 
gas (including the wide-spread, low-density and 
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sub-thermally excited component) in high-z galax- 
ies. The search for CO(l-O) in HzRGs at observ- 
ing frequencies below 50 GHz with the ATCA and 
EVLA will be a crucial complement to surveys of 
the higher CO transitions with ALMA. 

We thank Laura Pentericci for the use of her 
HST/NICMOS and VLA radio continuum images. 
Reproduced by permission of the AAS. The Aus- 
tralia Telescope is funded by the Commonwealth 
of Australia for operation as a National Facility 
managed by CSIRO. This research has made use of 
the NASA/IPAC Extragalactic Database (NED) 
which is operated by the Jet Propulsion Labo- 
ratory California Institute of Technology under 
contract with the National Aeronautics and Space 
Administration. 
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